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Abstract. We show on two examples, namely a calculation for charge symmetry breaking in 
pn —> dn that allows one to extract the quark mass difference induced part of the proton-neutron 
mass difference and a high precision calculation for pion-deuteron scattering and its implications 
for the value of the charged pion-nucleon coupling constant, how QCD tests can be performed from 
low energy hadronic observables. 
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INTRODUCTION 

In the last decades the effective field theory for low energy phenomena within the 
Standard Model, Chiral Perturbation Theory (ChPT), has developed to a mature tool 
to study hadronic phenomena at low energies with a clear cut connection to QCD — see 
Refs. [1, 2, 3] for recent reviews with emphasis on the %%, the single nucleon and the 
two-nucleon sector, respectively. 

One very useful application of ChPT is its use to extract from complex reactions more 
fundamental quantities that can be compared to QCD predictions straight forwardly. 
Those QCD predictions are calculated from first principles using lattice gauge theory 
techniques [4]. Since those are quite involved numerically the described interplay of 
effective field theory and numerical methods is very rewarding. For a long time hadronic 
reactions were studied using models. Although very successful in providing a qualitative 
picture of the reaction mechanisms, it is not possible to determine the accuracy of 
the calculation. Here effective field theories are in a clear advantage: since in their 
very nature they are controlled expansions in some small parameter, they allow for 
uncertainty estimates. This is why the role of model calculations is decreasing in recent 
years. The interplay of models, ChPT and lattice QCD is illustrated in Fig. 1. 



In this presentation two examples for the strategy outlined will be described. On 
the one hand a calculation will be sketched where from an analysis of the isospin 
violating forward-backward asymmetry of pn d7i + the quark mass induced piece 
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FIGURE 1. Illustration of different methods to compare properties of QCD to experiment. 



of the proton-neutron mass difference was extracted [6], and on the other hand a high 
precision calculation is discussed that allowed for an improved extraction of the 7iN 
scattering lengths [5] from an analysis of high accuracy pionic deuterium data. 



ISOSPIN VIOLATING NN -> d% 

Many of the isospin violating observables are dominated by effects from the pion mass 
difference, since, although m n o — m n + is typical for hadronic mass differences within 
isospin multiplets, isospin violating effects are enhanced due to the small pion mass. This 
is the reason for the importance of charge symmetry breaking (CSB) reactions (under 
charge symmetry up and down quark get interchanged), for here the pion mass difference 
does not contribute and therefore effects from different quark masses get enhanced. 

In this section we focus on the CSB null observable Aft (pn —> dn ) — the forward- 
backward asymmetry in pn — >» dn°. This is a null observable since in a charge symmetric 
world the final state fixes the total isospin to 1 and therefore the initial state gets projected 
on isospin 1 and is to behave as a proton-proton/neutron-neutron pair where forward and 
backward are not defined. At TRIUMF this observable was found to be [7] 

A fb = [17.2±8(stat.)±5.5(sys.)] x 1(T 4 . (1) 

In this section we will briefly describe the first complete NLO calculation for this re- 
action within ChPT. Also another CSB null-observable was measured recently, namely 
the total cross section for d d —> a7t° [8], however, since the four-nucleon dynamics in- 
volved in the reaction is a lot more complicated, no complete theoretical analysis exists 
yet for this reaction — see Refs. [9, 10] for some preliminary studies. 

The calculation described here became possible due to recent advances in developing 
a systematic power counting for reactions of the type NN —> NNk that was complicated 
by the presence of the large initial momentum p^ = y^m^M^, with m K (Af#) for the 



pion (nucleon) mass, which calls for a different expansion parameter, namely [11, 12] 




^prod — A 

where for the last identity the chiral symmetry breaking scale A x was identified with the 
nucleon mass. Nowadays the ChPT calculations have basically replaced the phenomeno- 
logical calculations (see Ref. [13] and references therein) that dominated the field before. 
For a recent review see Ref. [14]. The reason why Aft, is linked to the proton-neutron 
mass difference is that the transformation properties of the quark mass term in QCD 
under axial rotations dictates a link between mass differences of heavy hadrons and the 
isospin violating pion scattering off the very same hadrons [15, 16, 17]. For the case of 
the pn —> dn° this was first studied in Refs. [18, 19], however, these calculations were 
not complete, for besides diagram (a) of Fig. 2, where the mentioned isospin violating 
TlN scattering enters, also diagram (b) enters. It gives a non-vanishing contribution since 
the isospin conserving kN interaction is energy dependent and therefore gets sensitive 
to the different energy transfer in 7i + exchange, equal to m n /2 + M p —M n , and in 7i~ 
exchange, equal to m K /2 +M n — M p . In general CSB due to electromagnetism and due 
to quark mass differences enter with similar strength. Here, however, it so happens that 
the sum of diagram (a) and (b) are proportional only to 8M qm — the strong part of the 
proton-neutron mass difference. Using the results of Ref. [20] as input, the calculation 
revealed 

Al° = (11.5±3.5)xl0-*^. (2) 

The calculation sketched refers to a leading order calculation, however, all contributions 
at NLO, namely one loop diagrams with virtual photons, cancel [10] — the reason for 
this cancellation is now understood [21]. Thus the uncertainty was estimated to be of 
order Xp Y od ~ 30%. Using the experimental result of Eq. (1), Eq. (2) may be converted 
to give 

<5M qm = (1.5 ±0.8 (exp.)±0.5 (th.)) MeV , (3) 

where the first (second) uncertainty follows from the uncertainty of the experiment 
(calculation). In Fig. (3) this result is compared to previous extractions: one directly 
from the proton-neutron mass difference using the Cottingham sum rule to quantify the 
electromagnetic contribution [22] and one from lattice QCD [23]. Note, the calculation 
presented has a comparable accuracy to the other extractions — thus an improved 
measurement would be very desirable. 



ltd SCATTERING LENGTH AND ITS IMPLICATIONS FOR g nNN 

The problem of nd scattering has been studied theoretically already for many decades 
using phenomenological approaches, however, nowadays the high accuracy of modern 
experiments calls for improved tools for the analysis. Especially, a consistent treatment 
of strong and electromagnetic few-body effects is essential for a controlled extraction of 
the quite small isoscalar pion-nucleon scattering length a + , for especially electromag- 
netic effects might even outnumber its contribution to the nd scattering length [24]. 
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FIGURE 2. Leading order diagrams for the isospin violating s-wave amplitudes of pn — >> d7l°. Solid 
(dashed) lines denote nucleons (pions). Diagram (a) corresponds to isospin violation in the TtN scattering 
vertex explicitly whereas diagram (b) indicates an isospin-violating contribution due to the neutron- 
proton mass difference in conjunction with the time-dependent Weinberg-Tomozawa operator (see text 
for details). 
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FIGURE 3. Comparison of different extractions of the quark mass induced neutron-proton mass dif- 
ference. The points are from Refs. [22] (Direct extraction), [23] (Lattice), and [6] (this work). The inner 
(red) error bars on the last point refer to purely the theoretical uncertainty. 



There were various important advances that made a high accuracy calculation for the 
7td scattering length, reported in Ref. [25], possible: various subleading contributions 
were shown to vanish [26], there exists a calculation for nN scattering of the necessary 
accuracy [27], the role of various few body corrections is understood [28, 29, 30], the 
role of the nucleon recoils is understood [31, 32], and dispersive and Delta corrections 
are nowadays under control quantitatively [33, 34]. In this section we will briefly sketch 
the results of Ref. [25], with special emphasis on isospin violating parts. 

The data for hadronic scattering lengths is best deduced from high accuracy measure- 
ments of pionic atoms [35] together with properly improved Deser formulae [36, 37] — 
for a recent review see Ref. [38]. 

The theoretical limit for the accuracy of a calculation of the TlNN —> 7tNN transition 
operator is set by the first (AW) 2 ;r 2 -counter term. In the power counting of Ref. [30] it 
appears at &(% 2 ) relative to the leading two-nucleon operator shown in Fig. 4 (d\), with 
X = m K /M]sf. We thus aim at a calculation with up-to and including <^(# 3 / 2 ) — square 
root orders appear due to the connection between pion production (see previous section) 




FIGURE 4. Topologies for it d scattering. Solid, dashed, and wiggly lines denote nucleons, pions, and 
photons, respectively. The blobs indicate the deuteron wave functions. 



and the dispersive corrections [33] as well as the numerical proximity of the Delta- 
nucleon mass difference and introduced above [34]. The diagrams that contribute up 
to this order, besides those for Delta and dispersive terms, are shown in Fig. 4. Naively, 
one might expect the most important isospin violating contributions to 7id scattering 
from different pion masses in the leading contributions, especially in the diagrams shown 
in Fig. 4 (d\) and (dfo). However, it is the subtle interplay of one-nucleon and two- 
nucleon operators, driven by the Pauli principle, already discussed in Ref. [31], and the 
orthogonality of the nuclear wave functions [32], that strongly suppresses these effects. 

More difficult is the treatment of photon loops that might get enhanced due to the 
masslessness of the photon together with the smallness of £, the deuteron binding energy. 
For example, one finds for the leading contributions of diagrams ( d§) in an expansion in 
small momenta 

_ f JVW(P-qMp) 8tt a- ( ( fT 

4 ~ J q 2 (q 2 + 2m w (e + p 2 /M p ))^ 3^V^V W% 

and (ds) of Fig. 4 



_ r d^pd^HpMp) 8tt a- , 



where *F denotes the deuteron wave function. Individually the amplitudes appear to 
be enhanced by a factor ^Jrn n j£ ~ 8 compared to the dimension analysis estimate, 
however, in the sum the enhanced pieces cancel. Similar cancellations can be observed 
for the other potentially infrared enhanced contributions. 

At the end it turns out that most of the additional contributions cancel pairwise and 
thus already the leading diagram — Fig. 4 (d\) — largely exhausts the value of the 7id 
scattering length. The numerically most important corrections are provided by an isospin 
violating piece to the 7iN scattering length and the triple scattering diagram (ds) — from 
the dimensional analysis this diagram contributes at <^(# 2 ), it is, however, enhanced by 
a factor 7t 2 due to its special topology and thus needs to be considered [30]. In addition, 
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FIGURE 5. Combined constraints in the a + -a~ plane from data on the width and energy shift of nH, 
as well as the %D energy shift. 



from pionic atoms it is not possible to extract a + directly, but only the combination [24] 

1 



d + = a + + 



1 +m n /M N 
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with F K for the pion decay constant and the additional low energy constants c\ and f\ . 
The combined analysis for pionic hydrogen and pionic deuterium data yields from the 
la error ellipse (c.f. Fig. 5) 



a 



(1.9±0.8) • l(T 3 m-\ a~ = (86.1 ±0.9) • lO -3 ™" 1 , 



(4) 



with a correlation coefficient p a -a+ = — 0.21 . We find that the inclusion of the 7iD energy 
shift reduces the uncertainty of <3 + by more than a factor of 2. Note that in the case of 
the 7iH level shift the width of the band is dominated by the theoretical uncertainty in 
Aa n - p , whereas for the 7iH width the experimental error is about 50% larger than the 
theoretical one. For details on the error budget see Ref. [25]. 

Taken together with c\ = (— 1.0±0.3)GeV _1 [39] and the rough estimate \f\\ < 
1.4GeV _1 [40], Eq. (4) yields a non-zero a + at better than the 95 % confidence level: 



(7.6±3.1)-10" 



(5) 



The final result for a + is only a little larger than several of the contributions considered 
in our analysis. This emphasizes the importance of a systematic ordering scheme, and 
a careful treatment of isospin violation and three-body dynamics. A reduction of the 
theoretical uncertainty beyond that of the present analysis will be hard to achieve without 



additional QCD input that helps pin down the unknown contact-term contributions in 
both the 7iN and kNN sectors. 

As it was argued in the introduction, 7iN scattering lengths are interesting quantities 
by themselves, especially since they can be extracted from lattice QCD calculations 
relatively easily. In the last part of this section we will show that in addition they also 
provide an important link between pion-nucleon and nucleon-nucleon scattering and 
in this sense a non-trivial consistency check for our current understanding of these 
fundamental reactions. Extracted from a careful analysis of 7iN scattering data, for long 
the charged pion nucleon coupling constant was believed to be g 2 /4n = 14.2 ± 0.2 [41]. 
However, when extracted from NN scattering [42], the value deduced reads g 2 /47i — 
13.54 ±0.05, where the error includes only the fitting uncertainty and not any possible 
systematic uncertainties. The use of the work presented to resolve the question on the 
value of g c becomes explicit, when using the Goldberger-Miyazawa-Oehme (GMO) sum 
rule [43] 2 . It reads 

Here the 7tN scattering lengths, now known to higher accuracy (c.f. Fig. 5), appear as 
subtraction constants for the dispersion integral 

47T 2 J V^+^l 

that may be expressed in terms of observable cross sections. Values for this integral 
can be taken from Refs. [45, 46]. Combining the findings of these works gives for 
the integral (—1.073 ± 0.034) mb [39], which is consistent with previous extractions 
(see Ref. [47]). The GMO sum rule was used before to pin down the value of g c , 
however, different analyses came to different answers. While Ref. [45] found a value 
as large as g 2 /47i — 14.11 ±0.05 ±0.19, where the first uncertainty is statistical while 
the second is systematic, Ref. [47] found values for g 2 /47i between 13 and 13.3. Also 
other, more general analyses from this group reported lower values, namely g 2 /47i — 
13.75 ±0.15 [48], g 2 c /4n = 13.76 ±0.01 [49]. Here the uncertainties only represent 
the statistical uncertainty. No attempt was made to quantify the systematics. More 
recently Ref. [46] presented g 2 /47i — 13.56 ±0.36 from a GMO analysis. There were 
no further developments in the last three years and the key players basically 'agreed to 
disagree' [50]. The basic improvements provided by the analysis discussed are that for 
the first time isospin violating corrections were included completely and consistently 
and, as the result of using a systematic effective field theory, it became possible to 
properly control the uncertainties of the 7lN scattering lengths. With the in this way 
improved input we find 

g 2 c /4n = 13.69±0.12±0.15 = 13.7±0.2 , 
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2 In the definition of g c there appear subtle issues due to the Coulomb poles [44] as well as other infrared 
singularities. For a detailed analysis we refer to Ref. [39]. 



where the first error gives the uncertainty in the scattering lengths and the second in the 
integral. From our analysis we therefore conclude that a value for g*/An above 14 is 
largely excluded. 

SUMMARY 

Modern lattice QCD calculations allow for first principle calculations of hadronic ob- 
servables like scattering lengths [4]. However, in most cases those quantities are not 
directly accessible from experiments but need to be extracted from reactions with com- 
plicated few body dynamics. In this presentation on two examples, the extraction of 
the quark mass induced proton-neutron mass difference, AM qm , from the forward- 
backward asymmetry in pn —> d7t°, and the extraction of the pion-nucleon scattering 
lengths from data on pionic hydrogen and deuterium, it is demonstrated that ChPT can 
be employed to extract, with controlled uncertainty, the quantities of interest from com- 
plex reactions to allow for a comparison to lattice data and thus for non-trivial test of 
QCD dynamics at low energies. 
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